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Abstract 

Critical Point Symmetries (CPS) appear in regions of the nuclear chart where a rapid change 
from one symmetry to another is observed. The first CPSs, introduced by F. Iachello, were 
E(5), which corresponds to the transition from vibrational [U(5)] to 7-unstable [0(6)] be- 
haviour, and X(5), which represents the change from vibrational [U(5)] to prolate axially 
deformed [SU(3)] shapes. These CPSs have been obtained as special solutions of the Bohr 
collective Hamiltonian. More recent special solutions of the same Hamiltonian, to be de- 
scribed here, include Z(5) and Z(4), which correspond to maximally triaxial shapes (the 
latter with "frozen" 7 = 30°), as well as X(3), which corresponds to prolate shapes with 
"frozen" 7 = 0°. CPSs have the advantage of providing predictions which are parameter free 
(up to overall scale factors) and compare well to experiment. However, their mathematical 
structure [with the exception of E(5)] needs to be clarified. 
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1 Introduction 



Critical point symmetries [H Ej, describing nuclei at points of shape phase transitions 
between different limiting symmetries, have recently attracted considerable attention, since 
they lead to parameter independent (up to overall scale factors) predictions which are found 
to be in good agreement with experiment [3 EJ [7] . The E(5) critical point symmetry 
PQ is supposed to correspond to the transition from vibrational [U(5)] to 7-unstable [0(6)] 
nuclei, while the X(5) critical point symmetry 2\ is assumed to describe the transition 
from vibrational [U(5)] to prolate axially symmetric [SU(3)] nuclei. Both symmetries are 
obtained as special solutions of the Bohr Hamiltonian |HJ. In the E(5) case pQ the potential 
is supposed to depend only on the collective variable (3 and not on 7. Then exact separation 
of variables is achieved and the equation containing (3 can be solved exactly PU E] for an 
infinite square well potential in j3, the eigenf unctions being Bessel functions of the first 
kind, while the equation containing the angles has been solved a long time ago by Bes 
[TO] . In the X(5) case j2] the potential is supposed to be of the form u{0) + u{^). Then 
approximate separation of variables is achieved in the special case of 7 ~ 0, the /3-equation 
with an infinite square well potential leading to Bessel eigenfunctions, while the 7-equation 
with a harmonic oscillator potential having a minimum at 7 = leads to a two-dimensional 
harmonic oscillator with Laguerre eigenfunctions |2J. In both cases the full five variables 
of the Bohr Hamiltonian [H] (the collective variables (3 and 7, as well as the three Euler 
angles) are involved. The algebraic structure of E(5) is clear, since the Hamiltonian is the 
second order Casimir operator of E(5), which corresponds to the square of the momentum 
operator in five dimensions (see [TT|IT2*] and references therein), while an SO(5) subalgebra 
(generated by the angular momentum operators in five dimensions) exists. The algebraic 
structure of X(5) (if any, since X(5) is an approximate and not an exact solution) has not 
been identified yet. 

It is of interest to identify additional special cases leading to analytical solutions of the 
Bohr Hamiltonian, and to examine their relation to critical behaviour of nuclei, clarifying 
in parallel their algebraic structure. 

It has been known for a long time that the Bohr equation gets simplified in the special 
case of 7 = 30° |131 Hlj , since two of the principal moments of inertia become equal in this 
case, guaranteeing the existence of a good quantum number (the projection a of angular 
momentum on the body-fixed x' axis), although the nucleus possesses a triaxial shape. In 
other words, the Hamiltonian possesses a symmetry, while the shape of the nucleus does 
not. By allowing the potential to be of the form u((3) + 11(7), and by permitting 7 to vary 
only around 7 ~ 30°, approximate separation of variables is achieved [T^], similar in spirit 
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to the X(5) solution. The /3-equation with an infinite square well potential leads then to 
Bessel eigenfunctions, while the 7-equation with a harmonic oscillator potential having a 
minimum at 7 = 30° takes the form of a simple harmonic oscillator equation. The full five 
variables of the Bohr Hamiltonian are involved in this case, while the algebraic structure 
(if any, since the solution is approximate) is yet unknown. This solution, which has been 
called Z(5) [TH], is presented in Section 2. 

Separation of variables becomes exact by "freezing" the 7 variable to the special value 
of 7 = 30°, in the spirit of the Davydov and Chaban ^B] approach. Then the /3-equation 
with an infinite square well potential leads to Bessel eigenfunctions ^2]; while the equation 
involving the Euler angles and the parameter 7 (which is not a variable any more) leads to 
the solution obtained by Meyer-ter-Vehn ^3]. The projection a of angular momentum on 
the body-fixed x' axis is a good quantum number also in this case. Only four variables (0 
and the three Euler angles) are involved, while the full algebraic structure is yet unknown. 
It has been remarked ^2], however, that the ground state band of this model coincides with 
the ground state band of E(4), the Euclidean algebra in four dimensions. This solution, 
which has been labelled as Z(4) [T2j, is presented in Section 3. 

The question arises then of what happens in the case one "freezes" the 7 variable to 
the value 7 = 0, which corresponds to axially symmetric prolate shapes, for which the 
projection K of angular momentum on the body-fixed 5-axis is a good quantum number. 
It turns out that only three degrees of freedom are relevant in this case, since the 
nucleus is axially symmetric, so that two angles suffice for describing its direction in space, 
while the variable describes its shape. Separation of variables becomes exact [d, the 
equation with an infinite square well potential leading to Bessel eigenfunctions, while 
the equation involving the angles leads to the simple spherical harmonics. The algebraic 
structure of this model is yet unknown. This solution, which has been called X(3) J7j, is 
presented in Section 4. 

Finally, in Section 5 the present results are briefly discussed and plans for further work 
are exposed. 



2 The Z(5) model 



The original Bohr Hamiltonian jH] is 
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where (3 and 7 are the usual collective coordinates, while Qk (k = 1, 2, 3) are the components 
of angular momentum in the intrinsic coordinate system and B is the mass parameter. 

In the case in which the potential has a minimum around 7 = tt/6 one can write the 
last term of Eq. (JTJ) in the form 4(Qf + Q\ + Q\) — ?>Q\. Using this result in the Schrodinger 
equation corresponding to the Hamiltonian of Eq. (JTJ), introducing [T] reduced energies 
e = 2BE/h 2 and reduced potentials u = 2BV/h 2 , and assuming [2] that the reduced 
potential can be separated into two terms, one depending on (3 and the other depending on 
7, i.e. u((3,j) = u{(3) + u(j), the Schrodinger equation can be separated into two equations 



1 d ^ d + ' 

((3 2 ) sin 37 $7 c?7 



(2) 



(3) 



where L is the angular momentum quantum number, a is the projection of the angular 
momentum on the body- fixed x'-axis (a has to be an even integer [H]), (ft 2 ) is the average 
of (3 2 over £(/3), and € = 6/3 + e 7 . 

The total wave function should have the form 7, 6^) = ^L, a (/?)?7(7)Pf f where 
#i (2 = 1, 2, 3) are the Euler angles, V(Qi) denote Wigner functions of them, L is the 
angular momentum quantum number, while M and a are the eigenvalues of the projections 
of angular momentum on the laboratory fixed i-axis and the body-fixed f-axis respectively. 

Instead of the projection a of the angular momentum on the x'-axis, it is customary to 
introduce the wobbling quantum number [T^ ITS] n w = L — a. Inserting a = L — n w in Eq. 
(J2J) one obtains 



±±f± + 

p A d(3 H 8(3 4/3 



2 (L(L + 4) + 3n w (2L - n w )) + u(J3) 



where the wobbling quantum number n w labels a series of bands with L = 11 
4, . . . (with n w > 0) next to the ground state band (with n w = 0) [Hj . 
In the case in which u{(3) is an infinite well potential 



£,L,n w (P) 
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one can use the transformation [2] £(/3) = /3 3//2 £(/3), as well as the definitions [2] 6/3 
z = (3k p, in order to bring Eq. (@J) into the form of a Bessel equation 
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Figure 1: (a) Intraband and interband B(E2) transition rates in the Z(5) model ^H], nor- 
malized to the B(E2;2 1)0 — > 0i ]0 ) rate. Bands are labelled by (s,n w ), their levels being 
normalized to 2i o- The (2,0) band is shown both at the left and at the right end of the 
figure for drawing purposes, (b) Same for the Z(4) model [T2] . 

Then the boundary condition £(@w) — determines the spectrum 

£p;s,v — £/3;s,n w ,L = (k s , u ) 2 , k s ,v = X S ,v/ @W, (8) 

and the eigenfunctions 

M/3) = Zs >nv> M = Zs,<*M = c s , u f3- 3/2 J u (k s , u p), (9) 

where x S)U is the sth zero of the Bessel function J u (z), while the normalization constants 
c SjU are determined from the normalization condition f3 4 ^ v (f3)df3 = 1. The notation for 
the roots has been kept the same as in Ref. [2j, while for the energies the notation E S}JIwj l 
will be used. The ground state band corresponds to s — 1, n w — 0. We shall refer to 
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the model corresponding to this solution as Z(5) (which is not meant as a group label), in 
analogy to the E(5) P and X(5) [2] models. 

The 7-part of the spectrum is obtained from Eq. Q, which, in the case of a harmonic 
oscillator potential having a minimum at 7 = tt/6, takes the form of a simple harmonic oscil- 
lator equation ^] . Similar potentials and solutions in the 7- variable have been considered 
in jH]. The total energy in the case of the Z(5) model is then 

E{s, n w , L, rij) = E + A(x SjV ) 2 + Bn^, (10) 

where is the quantum number of the oscillator occuring in the 7-equation. 

The calculation of B(E2) transition rates has been described in detail in Ref. The 
resulting level scheme is shown in Fig. 1(a). Experimental manifestations of Z(5) seem to 
appear in the 192 ~ 196 Pt region 



3 The Z(4) model 
3.1 The Z(4) solution 

In the model of Davydov and Chaban ^H] it is assumed that the nucleus is rigid with 
respect to 7-vibrations. Then the Hamiltonian depends on four variables (/?, 6i) and has 
the form ^B] 



H 



2B 



/? 3 <9/T d(3 



Q 



1 3 
4^ ^ sin 2 (7 - 



In this Hamiltonian 7 is treated as a parameter and not as a variable. The kinetic energy 
term of Eq. (Jllj) is different from the one appearing in the E(5) and X(5) models, because 
of the different number of degrees of freedom treated in each case (four in the former case, 
five in the latter). 

Introducing PQ reduced energies e = (2B/h 2 )E and reduced potentials u = (2B/h 2 )U, 
and considering a wave function of the form ty({3,9i) = (f)((3)ip(9i), where 0+ ( i = 1, 2, 3) 
are the Euler angles, separation of variables leads to two equations 
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In the case of 7 = 7r/6, the last equation takes the form 



-(Ql + 4Ql + 4Ql)-\ 



This equation has been solved by Meyer-ter-Vehn 
with 



I 2L 
16tt 2 (1 + 5, 



a,0J 



m) = 0. 

the eigenfunctions being 

+ (-i)M^ 



(14) 



(15) 



A = A L ,« = L(L + l)--t^, (16) 

where // and a are the eigenvalues of the projections of angular momentum on the laboratory 
fixed i-axis and the body- fixed x'-axis respectively, a has to be an even integer [Tlj. As 
in the previous section, the wobbling quantum number, n w [TU HE] , is introduced at this 
point. 

The "radial" Eq. (fT2*j) is exactly soluble in the case of an infinite square well potential 
[Eq. (j5J]. Using the transformation 0(/3) = /3 _1 /(/3), Eq. (|T2j) becomes a Bessel equation 

4- I F 
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Then the boundary condition f(fiw) — 
(jHj), while the eigenfunctions are 

<K/3) = = ^,l(/3) 



determines the spectrum, which is given by Eq. 



(3 J v ( k sl> /3), 



2 



W j2 ( 



(19) 



where the normalization constant c is determined from the condition Jq w f3 3 <j) 2 (f3)d/3 = 1. 
The notation for the roots has been kept the same as in Ref. while for the energies the 
notation E s>niWjL will be used. The ground state band corresponds to s — 1, n w — 0. This 
model will be called the Z(4) model. 

The calculation of B(E2)s proceeds as described in Refs. [T21II3]. The resulting level 
scheme is shown in Fig. 1(b). Experimental manifestations of Z(4) seem to appear in the 
128 ~ 132 Xe region jT^j. One can easily see that the spectra of the ground state band and the 
(3\ band, as well as the related transitions, are very similar to the ones predicted by the 
E(5) model PUEh while for the 71 bands the odd levels are very similar, while the even 
levels exhibit opposite staggering |12j . A partial explanation of this behaviour is given in 
the next subsection. 
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3.2 Relation of the ground state band of Z(4) to E(4) 

The ground state band of the Z(4) model is related to the second order Casimir operator 
of E(4), the Euclidean algebra in four dimensions. In order to see this, one can consider 
in general the Euclidean algebra in n dimensions, E(n), which is the semidirect sum ^J] of 
the algebra T„ of translations in n dimensions, generated by the momenta Pj = —id/dxj, 
and the SO(n) algebra of rotations in n dimensions, generated by the angular momenta 

( d d \ 

symbolically written as E(n) = T n © s SO(n) [20] • The generators of E(n) satisfy the 
commutation relations 

[Pi, Pj] = 0, [Pi, L jk ] = i(5 tk P, - 5 l3 P k ), (21) 

[Lij, L k i] = i[8i k Lji + 5jiL ik — 5nLj k — 8j k Lu). (22) 

From these commutation relations one can see that the square of the total momentum, P 2 , 
is a second order Casimir operator of the algebra, while the eigenf unctions of this operator 
satisfy the equation 

1 1 d -r n - 1 l + Uj{UJ + :- 2) )F(r) = k 2 F(r), (23) 



in the left hand side of which the eigenvalues of the Casimir operator of SO(n), u(uj + n — 2) 
appear [21J. Putting F(r) = r^ 2 ~ n ^ 2 f(r) and u = u + (n — 2)/2, Eq. is brought into 
the form 

the eigenf unctions of which are the Bessel functions f(r) = J v (kr). The similarity between 
Eqs. (J2lj) and (HI) is clear. 

The ground state band of Z(4) is characterized by n w = 0, which means that a = L. 
Then Eq. (118)1 leads to v — L/2 + 1, while in the case of E(4) one has v = u + 1. Then 
the two results coincide for L = 2u, i.e. for even values of L. One can easily see that this 
coincidence occurs only in four dimensions. 

It should be emphasized, however, that neither the similarity of spectra and B(E2) 
values of Z(4) to these of the E(5) model, nor the coincidence of the ground state band 
of Z(4) to the spectrum of the Casimir operator of the Euclidean algebra E(4) clarify the 
algebraic structure of the Z(4) model, the symmetry algebra of which has to be constructed 
explicitly, starting from the fact that 7 is fixed to 30°, for which the Bohr Hamiltonian 
possesses "accidentally" a symmetry axis (the body-fixed x'-axis). 
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4 The X(3) model 



In the collective model of Bohr [S] the classical expression of the kinetic energy correspond- 
ing to (3 and 7 vibrations of the nuclear surface plus rotation of the nucleus has the form 

|H1 122 

1 3 R 

^ (25) 



z fc=i z 



where (3 and 7 are the usual collective variables, B is the mass parameter, 

J k = AB(3 2 sin 2 (7 -27^/3) 

are the three principal irrotational moments of inertia, and u' k (k = 1, 2, 3) are the com- 
ponents of the angular velocity on the body-fixed fc-axes, which can be expressed in terms 
of the time derivatives of the Euler angles <fi, 8, if) 



lu[ — — sin 9 cos if} <f> + sin if) 8, u)' 2 = sin 9 sin if) <fi + cos if> 9, u' 3 = cos 9<j> + if). (26) 

Assuming the nucleus to be 7-rigid (i.e. 7 = 0), as in the Davydov and Chaban approach 
[IS], and considering in particular the axially symmetric prolate case of 7 = 0, we see that 
the third irrotational moment of inertia ^3 vanishes, while the other two become equal 
(Ji = J2 = 3B(3 2 ), the kinetic energy of Eq. (|23j) reaching the form 



T 



i3^ 2 K + ^) + |/3 2 = | r 



3/3 2 (sin 2 #0 2 + 6) 2 )+/3 2 



(27) 



It is clear that in this case the motion is characterized by three degrees of freedom. In- 
troducing the generalized coordinates q± = 0, q<i = 9, and ^3 = /3, the kinetic energy 
becomes a quadratic form of the time derivatives of the generalized coordinates 
T = BJ2ij=i9ij QiQj/2, with the matrix having a diagonal form 



9ij 



f 3/3 2 sin 2 9 
3/3 2 
V 1 



(2f 



(In the case of the full Bohr Hamiltonian [Hj the square matrix g^ is 5-dimensional and 
non-diagonal [221 ES] ) 

Following the general procedure of quantization in curvilinear coordinates one obtains 
the Hamiltonian operator 
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Figure 2: Energy levels of the ground state [s = 1), (3\ (s = 2), and /5 2 (s = 3) bands of X(3) 
|17j . normalized to the energy of the lowest excited state, 2f, together with intraband B(E2) 
transition rates, normalized to the transition between the two lowest states, B(E2\ 2f — > 
0i). Interband transitions are listed in Table 1. 



where is the angular part of the Laplace operator 



1 d . n d Id 2 
sinf — + 



sm9 89 89 sin 2 9 8(j) 2 ' 
The Schrodinger equation can be solved by the factorization 

y(P,9^)=F((3)Y LM (9,<j ) ), 



(30) 



(31) 



where Ylm(9,4>) are the spherical harmonics. Then the angular part leads to the equation 



AoY/ 



L(L + 1)Y! 



LM\ 



(32) 



where L is the angular momentum quantum number, while for the radial part F{(3) one 
obtains 

j4^- L ^ + >-M m ^ (33) 
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Table 1: Interband B(E2\ Li — > Lf) transition rates for the X(3) model |17] . normalized to 

the one between the two lowest states, B(E2; 2f — > Of). 

U -> L/ X(3) Li -> L/ X(3) -L, -> L/ X(3) 

2 -> 2 : 164.0 

2 2 -> 4i 64.5 2 2 -> 2 : 12.4 2 2 -> 0i 0.54 

4 2 -> 6i 42.2 4 2 -> 4 X 8.6 4 2 -> 2 X 0.43 
6 2 ^8! 31.1 6 2 ^6i 6.7 6 2 -> 4 X 0.51 

8 2 -> 10! 24.4 8 2 -> 8i 5.5 8 2 -> 6 X 0.56 

10 2 -> 12i 19.9 10 2 -> 10i 4.7 10 2 -> 8 X 0.59 
3 ^2 2 209.1 

2 3 -> 4 2 92.0 2 3 -> 2 2 16.2 2 3 -> 2 0.67 

4 3 -> 6 2 65.3 4 3 -> 4 2 12.2 4 3 -> 2 2 0.47 
6 3 -> 8 2 50.9 6 3 -> 6 2 10.1 6 3 -> 4 2 0.52 

8 3 -> 10 2 41.6 8 3 -> 8 2 8.6 8 3 -> 6 2 0.57 

10 3 -> 12 2 35.0 10 3 -> 10 2 7.5 10 3 -> 8 2 0.61 



As in the case of X(5) [2], the potential in j3 is taken to be an infinite square well [Eq. (jSJ- 
In this case F(f3) is a solution of the equation 



d 2 2d /, 2 

+ "7T "TTT + fe 



L(L + 1) 

3/5 2 



F(/3) = 



(34) 



in the interval < (3 < (3w, where reduced energies e = k 2 = 2BE/K 2 j2j have been 
introduced, while it vanishes outside. 

Substituting F(f3) = (i~ x l 2 j{fl) one obtains the Bessel equation 



d 2 Id l2 v 2 
^7 + ^+ \k 2 - -= 



d/3 2 f3d(3 



ft 2 
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where v 



L(L+1) , 1 



+ 7, the boundary condition being jifiw) — 0. The solution of (p^jl . 



which is finite at /3 = 0, is then 



F sL (/3) = ^/?- 1 / 2 J,(^/?), 



(36) 



with fc S) „ = x SjU /Pw and e Si „ = fc 2 ^, where x S)U is the s-th zero of the Bessel function of the 
first kind J v {k s ^(3w) and the normalization constant c = (3^ J 2 +1 (x StV )/2 is obtained from 
the condition Jq F 2 l ((3) j3 2 df3 = 1. The corresponding spectrum is then 
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It should be noticed that in the X(5) case [,2j the same Eq. occurs, but with v = 



+ |, while in the E(3) Euclidean algebra in 3 dimensions, which is the semidirect 



sum of the T3 algebra of translations in 3 dimensions and the SO (3) algebra of rotations 
in 3 dimensions [201 > the eigenvalue equation of the square of the total momentum, which 
is a second-order Casimir operator of the algebra, also leads [HI [20! ^° E<1- (JSHJ), but with 



From the symmetry of the wave function of Eq. (|3ip with respect to the plane which 
is orthogonal to the symmetry axis of the nucleus and goes through its center, follows that 
the angular momentum L can take only even nonnegative values. Therefore no 7-bands 
appear in the model, as expected, since the 7 degree of freedom has been fixed to 7 = 0. 

B(E2) transition rates are calculated as described in Ref. [T7j. The resulting level 
scheme is shown in Fig. 2 and Table 1. Experimental manifestations of X(3) seem to occur 
in 172 0s and 186 Pt [Tj\. An unexpected observation ^Zj is that the (3\ bands of the N=90 
isotones 150 Nd, 152 Sm, 154 Gd, and 156 Dy, agree very well with the X(3) predictions. These 
N=90 isotones are considered to be very good examples of X(5) [SI [201 123 I2H 120] , but the 
spacing within their fa bands is about half of that predicted by X(5). 

5 Discussion 

It should be remarked that in all of the above mentioned cases the Bessel eigenfunctions 
obtained are of the form J u (kf3), with v being of the form v = A + (^7^) > where n is 
the number of dimensions entering in the problem, while A = L(L + l)/3 in the cases of 
X(3) HZj and X(5) R, A = [L(L + 4) + 3n w (2L - n w )}/A in the cases of Z(4) [12 and 
Z(5) [IS] , with n w = L — a being the wobbling quantum number [TH], and A = r(r + 3) in 
the case of E(5), with r being the seniority quantum number characterizing the irreducible 
representations of the SO(5) subalgebra of E(5) [1. In the corresponding ground state 
bands one has n w = and r = L/2. 

It should also be mentioned that all the /3-equations mentioned above are also soluble 
[3Til 1ST] if the infinite square well potential is substituted by a Davidson potential [32] of the 
form u(f3) = (3 2 + (3$ / (3 2 , where (3q is the minimum of the potential, the eigenfunctions being 
Laguerre polynomials instead of Bessel functions in this case. A variational procedure has 
been developed |33| 13*1] , in which the first derivative of various collective quantities is max- 
imized with respect to the parameter /3 , leading to the E(5), X(5), Z(5), and Z(4) results 
in the corresponding cases. The solutions corresponding to the Davidson potentials lead to 
monoparametric curves [OSj connecting various collective quantities, where agreement with 
experimental data is very good. 




u = L + \. 
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Concerning future work, the clarification of the algebraic structure of the exactly soluble 
models X(3) and Z(4), as a prelude for the understanding of the algebraic structure of 
the approximate solutions X(5) and Z(5), is a challenging problem. The construction of 
analytical models including the octupole degree of freedom [SB] and/or the dipole degree of 
freedom is also receiving attention. 
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